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O ,■ Abstract 

<D 

T I^ , The alternative version of Hamiltonian formalism for higher-derivative theo- 

ries is proposed. As compared with the standard Ostrogradski approach it has 
the following advantages: (i) the Lagrangian, when expressed in terms of new 
variables yields proper equations of motion; no additional Lagrange multipliers 
are necessary (ii) the Legendre transformation can be performed in a straight- 
forward way provided the Lagrangian is nonsingular in Ostrogradski sense. The 
generalizations to singular Lagrangians as well as field theory are presented. 



Introduction 

It is a long-standing problem whether and why it is sufficient to use in physics the 
Lagrangians containing only first order time derivative. This is the more intriguing that 
adding higher derivatives may improve our models in some respects, like ultraviolet 
behaviour [U [2] (in particular, making modified gravity renormalizable [3] or even 
asymptotically free [1]); also, higher-derivative Lagrangians appear to be a useful tool 
to describe some interesting models like relativistic particles with rigidity, curvature 
and torsion [5] Moreover, almost any effective theory obtained by integrating out some 
degrees of freedom (usually, but not always, those related to high energy excitations) of 
the underlying "microscopical" theory contains higher derivatives. One can argue that 
the effective theory, being an approximation to perfectly consistent quantum theory 
need not to be considered and quantized separately. However, we are never sure if our 
theory is the basic or effective one; therefore, it is important to know whether it is at 
all possible to quantize the effective theory in a way which would correctly reproduce 
some aspects of the microscopic one. 
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First step toward the quantum theory is to put its classical counterpart in Hamilto- 
nian form. Standard framework for dealing with higher-derivative theories on Hamilto- 
nian level is provided by Ostrogradski formalism [6]- [10]. The main disadvantage of the 
latter is that the Hamiltonian, being linear function of some momenta, is necessarily 
unbounded from below. In general, this cannot be cured by trying to devise an alter- 
native canonical formalism. In fact, any Hamiltonian is an integral of motion while it 
is by far not obvious that a generic system described by higher-derivative Lagrangians 
posses globally defined integrals of motion, except the one related to time transla- 
tion invariance. Moreover, the instability of Ostrogradski Hamiltonian is not related 
to finite domains in phase space which implies that it survives standard quantization 
procedure (i.e. cannot be cured by uncertainty principle). 

Ostrogradski approach has also some other disadvantages. There is no straightfor- 
ward transition from the Lagrangian to the Hamiltonian formalism. In fact, Ostro- 
gradski approach is based on the idea that the consecutive time derivatives of initial 
coordinate(s) form new coordinates (& ~ q^~ l \ It appears then that the Lagrangian 
cannot be viewed as a function on the tangent bundle to coordinate manifold because 
it leads to incorrect equations of motion. Also, the Legendre transformation to the 
cotangent bundle (phase space) cannot be performed. One deals with this problem 
by adding Lagrange multipliers enforcing the proper relation between new coordinates 
and time derivatives of the original ones. This results in further enlarging of coordi- 
nate manifold; moreover, the theory becomes constrained (in spite of the fact that the 
initial theory may be nonsingular in the Ostrogradski sense, c.f. eq. (11.2J) below) and 
the Hamiltonian formalism is obtained by applying Dirac constraint theory, i.e., by 
reduction of the cotangent bundle to submanifold endowed with sympletic structure 
defined by Dirac brackets. 

In the present paper an alternative approach is proposed. It leads directly to the 
Lagrangians which, being a function on the tangent manifold, gives correct equations of 
motion; no new coordinate variables need to be added. Furthermore, for Lagrangians 
nonsingular in Ostrogradski sense the Legendre transformation takes the standard form. 
Our approach is also applicable to the most interesting case of singular Lagrangians 
(for example, those defining f(R) gravities [TT]). 

The paper is organized as follows. In Section [1] we consider nonsingular Lagrangians 
containing second and third order time derivatives. Constrained theories are discussed 
in Section [21 The general formalism is applied to mini-superspace formulation of f(R) 
gravity [12] in Section [3J In Section HJ the modifications necessary to cover the field- 
theoretic case are given. In Appendix we describe (for one degree of freedom) the 
generalization of our formalism to Lagrangians containing arbitrary high derivatives. 

1 Nonsingular Lagrangians of second and third or- 
der 

In this section we consider the Lagrangians containing second and third time deriva- 
tives which are nonsingular in Ostrogradski sense. Ostrogradski approach is based 



on the idea that the consecutive time derivatives of the initial coordinate form a new 
coordinates, g« ~ q^~ l \ However, it has been suggested [13]-[T6], [21] that one can 
use every second derivative as a new variable, qi ~ q i - 2% ~ 2 \ We generalize this idea 
by introducing new coordinates as some functions of the initial ones and their time 
derivatives. Our paper is inspired by the results obtainded in Ref. 



1.1 The case of second derivatives 

Let us start with Lagrangians containing time derivatives up to the second order, 

L = L(q,q,-q); (1.1) 

here q = (q 11 ), fi = 1, . . . , N denotes the set of generalized coordinates. The nonsingu- 
larity condition of Ostrogradski reads 

detf-?-^)/0. (1.2) 

\dq^dq v J 

In order to put our theory in the first-order form we define new coordinates gf , q 2 : 

q" = ql r = q?, q" = X^Qu <7i, <fe), (1-3) 

where x M are the functions specified below. 
We select an arbitrary function 

F = F(q 1 ,q ll q 2 ) 1 (1.4) 

subjected to the single condition 

Hwh)* a - (15) 

Now, x* 1 are defined as a unique (at least locally due to (11.21) ) solution to the following 
set of equations 

dL(q u qi,X) , dF(q 1 ,q 1 ,q 2 ) 

w + — M — = - (L6) 

The new Lagrangian, which is now a standard Lagrangian of first order, is given by 



dF(q 1 ,q 1 ,q 2 ) .^ 

dqf 
9F( q 1 ,q 1 ,q 2 ) . M dF(q 1 ,q 1 , q 2 ) 
dq 2 2 dq£ 



£{qi,qi,q2,q2) = L(q\,qi,x{qi,qi,q2)) + *V 2 4i 

^ g 2 + ^u x {qi, Qi,q2)- U-O 



It differs from the initial one by an expression which becomes "on-shell" a total time 
derivative. 



The equation of motion for q% yield 

d 2 F 



(X u -Qi) = 0, (U 



d<iidq% 
which, by virtue of (jl.5p . implies 

9 M = x M (9i, 91,92). (1.9) 

For the remaining variables q± one obtains 

dL d ( dL\ d 2 ( dL , 

0, (1.10) 



dqi dt \d<ji J dt 2 \dx 

and taking into account (jl.9p one gets the initial Euler-Lagrange equations. 

It is worth to notice that, contrary to the original Ostrogradski approach, the 
formalism presented above leads directly to the standard picture of Lagrangian as a 
function defined on the tangent bundle to coordinate space (with no need of enlarging 
of the latter by adding the appropriate Lagrange multipliers). 

Our Lagrangian (jl.7p is nonsingular in the usual sense so one can directly pass to 
the Hamiltonian picture by performing Legendre transformation leading to canonical 
dynamics on cotangent bundle. 

To this end we define the canonical momenta 

_ dC dL d 2 F .„ 8 2 F „ 8 2 F .„ 8F 
M dqi dtji dq^dqi 1 dcftdq" dq^dq^ 2 dq%' 

_ dC dF(q 1 ,q 1 ,q 2 ) 

P2^ = 1 rrn = 5-5 • (1.12) 



dq% dq: 



By virtue of (jl.5p the second set of equations can be uniquely solved (at least locally) 
for $ 

Q.i =qi(qi,Q2,P2)- (1-13) 

As for the first set (II. lip , we note that q^ appears (linearly) only in the fourth term 
on the RHS. Again, the same condition (jl.5p allows us to solve fll.lip for q^ , 

0.2 =92(91, 92, Pi, P2). (1-14) 

The Hamiltonian H is computed in standard way and the final result reads 

dF dF 

H = P^-L- w ti- W X^ (1.15) 

where everything is expressed in terms of qi, q 2 ,Pi and p 2 - We have checked, by direct 
calculation, that the canonical equations following from H are equivalent to the initial 
Lagrangians ones. 



There exists canonical transformation which relates our Hamiltonian to the Ostro- 
gradski one. It reads 

9i = Ox, 

0.2 = f l (qi,Q2,P2), . . 

V2p, = —§jji(qi,f,q2), 

where tilde refers to Ostrogradski variables and f^iqi, 92,^2) solve eqs. (II. 121) . i.e. 
f 1 = 0i{0ii02iP2)- The corresponding generating function $ has the form 



®(Qi,Pi, to,®) = <^Pin + F(,Qi,Q2,Q2)- (1-17) 

However, it should be stressed that Ostrogradski Hamiltonian is singular in the sense 
that the inverse Legendre transformation cannot be performed (contrary to our case). 
This means that the structure of sympletic manifold (phase space) as a cotangent 
bundle to coordinate manifold is not transparent if Ostrogradski variables are used. 
Let us conclude this part with a very simple example. The Lagrangian 

L = \e ltt ,Fq" + -(??, /^0, fi,v= 1,2 (1.18) 

is nonsingular in Ostrogradski sense provided (3 7^ 0. We take 

F = aqfqo 1 , a^O. (1.19) 

Then 

X» = ~^ v q\--^ (1.20) 

and 

C = -^(qtf ~ ^) 2 ~ y^« + ««• (1-21) 

Finally, the Hamiltonian reads 

1 A 2 A a 2 

h = -v^V2v + ^q(p^) 2 + -n e ^P2^q u 2 + ^(^) 2 - (1-22) 

It depends on an arbitrary parameter a. One can pose the question whether any 
relevant physical quantity may depend on a. The answer is no: all physical quantities 
are a- independent. Formally this can be shown using eqs f 1 1 . 1 6 j) and (11.171) . Indeed, 
the function generating the canonical transformation to Ostrogradski variables reads 

^{OliPl^O^OT) = OlPlp + &0.2 0.2- ( L23 ) 

The corresponding canonical transformation takes the form 

Pin =Plp, 

S-% (L24) 

02 - aP2pi 

P2n = aq.2 ', 



when inserted into the Hamiltonian (jl.22p it yields the standard Ostrogradski Hamil- 
tonian 

H = Pi^t + ^fe,) 2 - ~^p 2v + ^m 2 - (1-25) 

It does not depend on a. Therefore, the energy (energy spectrum in quantum the- 
ory) does not depend on a. The role of our a-dependent modification is to provide the 
formalism which yields standard Lagrangian dynamics and regular Legendre transfor- 
mation. 

The above explanation is slightly formal. We shall now look at the problem of 
a dependence from a slightly different point of view. Let us note that the classical 
state our system is uniquely determined once the values of q(t), q(t), q(t), q(t) at some 
moment t are given. Moreover, most physically relevant quantities are constructed via 
Noether procedure (they are either conserved or partially conserved, i.e., their time 
derivatives are defined by transformation properties of symmetry breaking terms in 
the action). As such they are expressible in terms of q, q, q and q. Therefore the latter 
are the basic variables. One can find their quantum counterparts provided we compute 
the relevant Poisson brackets. 

To this end we write out the canonical equations of motion following from eq. (jl.22p : 

•A 4 1 

Pl/j, = 0) 

_ A a 2 A* 

Pip, — p e (J.vP2v £f<?2 ■ 

They lead to the following relations 

(f = tf, 

•^ _ A , „ a u [ i -^'> 

q - ^nvP2v - -p(l2 > ' 

"n^ 2A 2 „ 2Aa v 1 

One can now find the Poisson brackets among q, q, q and q. The nonvanishing ones 
read 

{y , q} = -^<W {<f> q) = <W 

{q^q} = -p^ {$",«"} = f^, (1.28) 

{q",q U }=P^ {q":q} = ^e, v . 

Note that they are a- independent. Upon quantizing we get four observables obeying 
a-independent algebra. Any other observable including energy can be constructed out 
of them so its spectrum and other properties do not depend on a. 

1.2 The case of third derivatives 

Let us consider a nonsingular Lagrangian of the form 

L=L{q,q,q,q). (1.29) 
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It is slightly surprising that this case (and, in general, the case when the highest time 
derivatives are of odd order - see Appendix) is simpler. We define the new variables 

q» = ql ? = $, ? = (£, r = q%. (1.30) 

Next, the function F(qi,q 1 ,q 2 ,q 3 ) is selected which obeys 

det G£)^ 0; (L31) 

here q 3 are additional variables. The modified Lagrangian reads 

r( . . ■ s ( . . . dF(q 1 ,q 1 ,q 2 ,q 3 ) dF(q 1 ,q 1 ,q 2 ,q 3 ) IjU 
C{qi,qi,q 2 ,q2,q3,q3) = L(q l ,q 1 ,q 2 ,q 2 ) + -^ <fi + -^ q% 

dF(q 1 ,q 1 ,q 2 ,q 3 ) .^ dF(q 1 ,q 1 ,q 2 ,q 3 ) M 

It can be easily shown that the Euler-Lagrange equations for L yield the initial equa- 
tions for the original variable g M = q±. Again, as in the second-order case, the Legendre 
transformation can be directly performed due to the condition (11.311) . The momenta 
read 

dL d 2 F .„ d 2 F v d 2 F .„ d 2 F .„ OF 
dL dF 

P ^ = W, + H' (L34) 



dF(q 1 ,q 1 ,q 2 ,q 3 ) 

dq 3 



P^ = ^ • ( L35 ) 



By virtue of ( 11.311) one can solve ( 11. 35ft for q^, 

4*1= Qi(Qi,Q2, q3,Ps)- (1-36) 

Inserting this solution into eq. (II. 34ft one computes 

q.2 =q2(qiiq2,q3,P2,P3); (1.37) 

the solution is (at least locally) unique because L is, by assumption, nonsingular in 
Ostrogradski sense. Similarly, (II. 33ft can be solved in terms of q 3 : 

q3=qs(qi,q2,q3,Pi,P2,P3)- (1.38) 

Finally, Hamiltonian is of the form 

OF dF dF 

H = Pmtf + V2A2 -L- -j^tf - T^q^ - ^2 , (1.39) 



where everything is expressed in terms of q's and jo's (the terms containing q 3 cancel). 
As above, we have checked that the canonical equations of motion yield the initial equa- 
tion. The canonical transformation which relates our formalism to the Ostrogradski 
one reads 



~u ri 

^ = ^ ' 

<?3 = 0.2 > 

Pip =PiiM-^(Qi,f(Qi,Q2,Q3,Pa),Q2,Q3), V" 40 ' 

f>2p, = —§F(quf(qi,<k,q3,Pa),<b,qa), 

Ps M = P2p - §^(qi, f(qi, q2, q3,Ps), Q2, qs), 

where f 1 is the solution of eq. ( 11.351) . i.e., / M = q%. The relevant generating function 
reads 

®fauPi,Qi,ch,q3,ihi) = qiPiix + q2P3p + F(q 1 ,q 2 ,q 2 ,q 3 ). (1.41) 

Again, the advantage of our Hamiltonian over the Ostrogradski one is that the former 
is nonsingular in the sense that the inverse Legendre transformation can be performed 
directly. 

1.3 The second order Lagrangian once more 

By comparing Section 11.11 and 11.21 we see that the modified Hamiltonian formalism is 
somewhat simpler in the case of third order Lagrangian (actually, as it is shown in 
Appendix, this is the case for all Lagrangians of odd order). Namely, in latter case 
no counterpart of the condition (\l.Q\i is necessary. This will appear to play the crucial 
role in the case of singular (in Ostrogradski sense) Lagrangians (see Section [2] below) . 
Therefore, as a preliminary step, we consider here the second order Lagrangians as 
a special, singular case of third order ones. The resulting Hamiltonian formalism is 
then constrained. However, with an additional assumption that the function F does 
not depend on q% , one can perform complete reduction of phase space obtaining the 
structure described in Section [TTT1 
Let 

L = L(q,q,q), (1.42) 

and F = F(qi,qi,qs) obeys (11.3ip . We define 

n ■ ■ \ ri ■ \ j_ dF (quQi,q3) -p 

£(31,91,92,93,93) = M9i, 9i,92) + ^ 9i 

09" 



dF( q 1 ,q 1 ,q 3 ) mfi dF(q 1 ,q 1 ,q 3 ) ^ 
dq» %+ d£ 



^ \H\-,HiiH-i .a . u± vyi,yi,y3; » /-, AO , 

+ JO «3+ 7^ &■ ( L43 ) 



The relevant momenta read 



dL_ d 2 F .„ d 2 F „ d 2 F .„ dF_ 
VXil ~ dtf + dq'dqf 1 + d(f x dqf 2 + dtfdqf 3 + d<£> { ' 

P2» = 0, (1.45) 

OF 

n, = w - (1.46) 



There is one set of primary constraints (11.451) . On the other hand, due to the condition 

(jl.3ip gf and q 3 can be expressed in terms of qi,q2,q3,Pi,Ps- The Dirac Hamiltonian 

takes the form 

dF .„ _ dF 

dtf qi ~dtf 

where c M are Lagrange multipliers enforcing the constraints $1^ = p 2fl ~ 0. 
The stability of primary constraints implies 



H = Pi»q? -L- j-ti - —^ + c»p 2 „ (1.47) 



dL(q 1 ,q 1 (q 1 ,q 3 ,p 3 ),q 2 ) dF(q 1 ,q 1 (q 1 ,q 3 ,p 3 ),q 3 ) 
dq$ + dtf 



q ^ A _ fo __ ^Wb'/iWi^i^i^; ^ wi^iWii^iWi^; /-, 4g x 



In order to check the stability of secondary constraints $> 2fl we note that, as it can be 
verified by direct computation, 

{#,«} = 0. (1.49) 

Using (11.491) together with 

0^$ 2Ai = {$ 2M ,#}, (1.50) 

we arrive at the following condition 

d 2 L v d 2 L .„ d 2 L v 0L_ 5F_ n 

dq^dqf + dqtdq^ 1 + dq^dqf 2 + ' ^ fl# 3g? ~ U ' ( & j 

The initial Lagrangian is nonsingular and eq. fll.51 j) can be used to determine the 
Lagrange multipliers c v uniquely. Therefore, the are no further constraints. 

In order to convert our constraints into strong equations we define Dirac brackets. 
To this end we compute 

{&„,&„} =0, (1-52) 

<*»*"> = -agfe = - w ''" (L53) 



Moreover, 




f dL 


dL 


\dq$- 


dql 


which implies 





dF dF) „ f OF dL) d 2 L 



dql ' dq\ J \ <9gf ' dq 2 J dq 2 dcft ' 



;i.54) 



{ ^ m = 4m-^m =v " v - (L55) 

By assumption, W is a nonsingular matrix. Consequently, 



is also nonsingular and 



c _, = ( ™- y , _ (L57) 
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Dirac bracket takes the following form 

{;-}d = {v}-{-,U^~ 1 W 1 W^.-} 

-{•, ^}(rV{^, •} + {■, 02 M }(W- 1 ) AB/ {0i I/J •}. (1.58) 

The constraints $1^ depend on p 2fM only. We conclude from (I1.58J) that the Dirac 
brackets for q%, q^Pi^^p^ take the canonical form. Moreover, p 2ii = while q 2 can be 
determined from (jl.48p . Note that the solution for q 2 , by virtue of eq. fll.6]) reads 

(& = X ti (qi,4i(Qi,<l3,P3),q3)- (1-59) 

So, up to renumbering q 2 -H- ^3 we arrived at the same scheme as in Section 11.11 
In order to illustrate the above approach, we use the same example as before: 

L = Xe^q + ^)\ P^O, (1.60) 

and 

F = aq?q%, a^O. (1.61) 

Then H takes the form 

H = -Pi m P3 M e^paptg - ~(q 2 ) 2 ~ aq^q^ + —e^ v p 2tl q 2 - -pi^jty, (1.62) 

a a I p p 



while the constraints are 



h(i = P2 M , 4>2f, = e^ u p 3l/ + j3q% + aq%, (1.63) 

a 



/" 



and serve to eliminate p 2il and q 2 

P2^ = 0, q% = —e^psv - -qg. (1.64) 

ap p 

Inserting this back into the Hamiltonian we arrive at the following expression 

H = -Pitfto? + 2tfQ^) 2 + -g e ^P^3 + Js^) 2 ( L65 ) 

which coincides with the one given by eq. (11.22)) provided the replacement q 2 -H- q 3 , 
P2 -H- P3 has been made. 

2 Singular Lagrangians of the second order 

In this section we consider the second order Lagrangians 

L = L(q,q,q), (2.1) 
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which are singular in the Ostrogradski sense, i.e. 

det(W%,) = det (^f) = 0. (2.2) 

For standard Ostrogradski approach to such singular Lagrangians see, for example, 

Ref. [H1I2S|. 

The formalism of Section 11.11 is not directly applicable because due to eq. ( 12.21) , 
eqs. (11.61) cannot be solved to determine the functions x^ ■ Moreover, eqs. (11.61) put in 
this case further restrictions on the form of F. 

In order to get rid of these problems we will follow the method of Section 11.21 and 
consider L as a third order singular Lagrangian. From this point of view its singularity 
comes both from eq. (12. 2p and the fact that the third order time derivatives are absent. 
Given a singular Lagrangian L we select a function F = F(q 1 ,qi,q 3 ) obeying (ll.3ip 
and define 

M9i, 91,92,93,9s) = M9i, 9i,92) + ^ 9i 

0<R 



dF( q 1 ,q 1 ,q 3 ) .^ dF(q 1 ,q 1 , q 3 ) ^ 
dq» Q3 + dqt 



aJi ? 3 H -^u 92- {*•*) 



As before, the canonical momenta given by (11.451) provide the primary constraints 
while (11.441) and (11. 46ft allow us to compute eft and q 3 . The Hamiltonian is given by 
eq. A1.47J1 . The secondary constraints read again 



^(9i, 91 (91,93,^3), 92) ^(9i, 9i(9i,9 3,P3), 93) 
dq 2 d Qi 



n^^s _ ^-^wi,yHyi,y3,^3;,y2; . ^ iyi,yuyi,y3,K3;,y3; ,« ,s 

u ~ ^ 2 m 53* + a^* ' V A > 



Now we have to investigate the stability of $2^- To this end we assume that W has 
rank K, K < N; this implies the existence of J = N — K linearly independent null 
eigenvectors 7^(91, 9i, 92), o = 1, 2, . . . , J, 

W, ul l = 0. (2.5) 

Equations ( I1.51J) do not determine uniquely the Lagrange multipliers c M ; on the con- 
trary, we get new constraints of the form 

° = ^ = 7 « [age** + Ww q * + n » ~ H ~ H ) ; (2 ' 6) 

here, as previously, eft = gf(gi, 93,^3), so the above constraints contain qi, q 2 , 93,^1 and 

We have started with third order formalism; therefore, our phase space is 6N- 
dimensional. As in nonsingular case (Section [T|) we would like to eliminate q 2 s and 
p 2 s. The latter are equal zero by primary constraints <3>i M . As far as q 2 s are considered 
the situation is more involved. 
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First, by virtue of the assumption (j2.2p about W we can determine from eqs. (12.41) 
K variables q 2 in terms of qi,Pi,q3,Pa and the remaining q 2 s. By substituting the 
resulting expression back to eqs. (12 Ah we arrive at J constraints on qi,pi,q 3 and p 3 . 
We denote these new constraints by tfj a (qi,q 3 ,pi,p 3 ). Let us now concentrate on the 
constraints (12. 6p . In general, they contain the q 2 variables and imply the constraints 
on qx,q3,Px,P3 only provided q 2 s enter in the combinations which can be determined 
from eqs. (j2.4p . In order to decide if it happens consider the variations 5q 2 which do 
not change the RHS of ( 12.41) . From the definition of W M „ we conclude that such 8q 2 
are linear combinations of 7^ (see (12.51) ) . If the RHS of ( 12.61) are stationary under such 
variations 5q 2 , eqs. (12.41) and ( 12.61) can be combined to yield the constraints which do 
not depend on q 2 s. The relevant condition reads 

3o - 7 f = 3 6=1,..., J; (2.7) 



A< 



dq 

where a takes M values, which without loss of generality can be chosen as a = 1, ... , M. 
In this way we obtain M new constraints on qi,p±, ^3,^3- 
One can check that 

^ „ = „ „ ( d 2 L d 2 L \ 

dq» lb lhla \dqtdql dftdtf J " l ^ 

By virtue of (11.551) we find 

Wa, M « 7^{02 M , K}, a = 1, . . . , M, b = 1, . . . , J. (2.9) 

Let us summarize. For the nonsingular second order Lagrangian viewed as a singular 
third order one, (qi,pi, 93,^3) forms the reduced phase space; no further constraints 
exist. On the contrary, in the singular case qi,pi, qz,Pz are still constrained. First, there 
exist J constraints '0a(<Zi>Pi> 93, P3) j moreover, if some (say - M) ^'s are in involution 
(on the constraint surface) with all ip's there exist additional M constraints following 
from eqs. (I2.4[) and (12.61) . This agrees with the conclusions of Ref. [17] . 

In general, for singular Lagrangian it is not possible to determine uniquely all 
Lagrange multipliers & '. However, we are in fact interested only in dynamical equations 
for qi, q^,pi and p^. Therefore, we can use the following Hamiltonian 

dF OF 



On the constraint surface it does not depend on q 2 , s, 

OH dL dF 



dq^ dq^ dqt 



0. (2.11; 



The existence of further secondary constraints depend on the particular form of the 
Lagrangian. 

Finally, let us note that the canonical transformation ( ll.40p leads to the form 
of dynamics presented in Ref. [T7J. However, within our procedure the Legendre 
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transformation from the tangent bundle of configuration manifold to phase manifold 
is again straightforward (if one takes into account standard modifications due to the 
existence of constraints). 

Singular higher derivative Lagrangians were also considered in [18]. The authors 
considered the physically important case of reparametrization invariant theories (higher- 
derivative reparametrization invariant Lagrangians appear, for example, in the descrip- 
tion of radiation reaction [19]). In their geometrical approach the image of the Legendre 
transformation form a submanifold of some cotangent bundle. This suggests that in 
the case of higher-derivative singular theories it is advantageous to start with enlarged 
phase space; this agrees with our conclusions. 

To conclude this section with a simple example consider the following Lagrangian 

L = Xe^q u + ^(q i ) 2 , fi,v=l,2. (2.12) 

It is singular and the matrix W (eq. (I1.53P ) is of rank 1 for (3 ^ and for (3 = 0. We 
take F as 

f = a«. (2-13) 

Assume first j3 ^ 0. Then 

C = Xe^qfq" + § fe 1 ) 2 + «« + «M, (2-14) 

and 

Pin = XepvqZ + aqg, 

P2» = 0, (2.15) 



P^ = aq_x- 



The primary constraints are 



while the Hamiltonian reads 



$1^=^-0; (2.16) 



H = -pinPSfi e»vP3»q% ~ o {(ilf ~ «?2^3 + c il P2 li - (2.17) 

a a I 



One easily derives the secondary constraints 

« $ 21 = ±p 32 _ p q \ - aql 

« $22 = ±P3i + aql { ' ' 

The stability for $ 2m yields 

« {$21, #} = 2Xqj -/3c 1 -p u , (2.19) 

« {$ 22 , if} = p 12 + 2Ag 2 1 = $ 3 . (2.20) 

Equation ( 12.191) allows us to compute c 1 , 

c 1 = i(2Ag 2 2 -p 11 ), (2.21) 
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while (I2.20p provides a new constraint. Its stability enforces c 1 = which together 
with ( 12.21 j) yields further constraint 

0^$ 4 = ^(2Ag 2 2 - Pll )- (2.22) 

Finally, differentiating the above equation with respect to time we get c 2 = 0. The 
resulting Hamiltonian is 

H = -P1^3 M <>P3/^2 - o (?2) 2 - «?2^3- ( 2 - 23 ) 

a a I 

Still we have to take into account the constraints $2 M , ^3 an d $4- The latter two can 
be rewritten as 

$3 M = Pi? - 2Xe^ql (2.24) 

$2^ and $3^ are now used in order to eliminate all variables except q^pi^ and q^,P32- 
The only nonstandard Dirac bracket reads 

{dlpm}D = \. (2.25) 



The Hamiltonian, when expressed in terms of unconstrained variables, takes the form 

H = -P12P32 ~ -rPnqj + ^7^(Pi2) 2 - (2.26) 



Let us note that the limit /3 — > is smooth. Of course, we could put (3 = from the 
very beginning and arrive at the same conclusion. 

3 An example: mini-superspace formulation of f(R) 
gravity 

As a more elaborate but still a toy example we consider mini-superspace Hamiltonian 
formulation of f(R) gravity [12]. We consider the following (LFRW - type) metrics 

ds 2 = -N 2 dt 2 + a 2 dx 2 . (3.1) 

Under such reduction the Lagrangian of f(R) gravity takes the form 

L(a,N) = ±Na 3 f{R), (3.2) 

where the curvature is given by 
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We see that L depends on second time derivatives. We proceed along the lines described 
in Section [TJ The basic dynamical variables are chosen as follows 

ax = a , di = a, Ni = N, N t = N, a 2 = R, (3.4) 

while 

a = x(ai, di,iVi,iVi,02), (3.5) 

is determined by eq. ( 11. 6ft once appropriate F is selected. We take 

F = -3o?/'(o 2 )di; (3.6) 

under the assumption /" ^ 0, eqs. (13.61) and (11.61) yield 

a 2 = R. (3.7) 

Solving ( 13. 7p with respect to a we find 



R-— T ^((2-N 1 )al-a 1 a l N 1 )). (3. 



6 V N Wx 
The modified Lagrangian reads 

1 ■**. n/ n w/ n / 9 6aia? 1 ,„ T 3a?diiV 



— 3/"(a2)a^did2. (3.9) 

It is straightforward to check that £ leads to the correct equations of motion. In order 
to simplify our considerations we introduce new variable 

711 = ^1/(02). (3.10) 

In terms of new variable £ reads 

r * 3 /O2) n . 2 -// x , 6oidf/ 2 (a 2 ) 

2 f'(a 2 ) m 

1 3 3a?didi/(a 2 ) , . 
- -o?nia 2 . (3.11) 

2 Til 

Now, we compute the canonical momenta: 

Pi = ^ = -^/(oaK (3.12) 

OTl! Til 

5£ is • *//- ^ , 12 aidi/' 2 (a 2 ) Q a 2 ,,, v , qiq , 

7Ti = —- = -I801O1/ (a 2 ) H 3—/ (a 2 )n 1; (3.13) 

aai ni n\ 

*=fj = 0. (3.14, 
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One can solve f l3.12[) and (13. 13[) in terms of di and hi. We form the Hamiltonian 

rhPiTh _ n 1 a\f(a 2 ) n\p\ _ 2n x p\ 

3a 2 J'(a 2 ) 2f'(a 2 ) a 3 J'(a 2 ) 3a? 

+ -ala 2 rii + pn 2 = H + /nr 2 . (3.15) 

Now, we investigate the stability of $1 = 7r 2 constraint 

The stability condition for $ 2 determines /i; an explicit expression for [i is irrelevant 
for what follows. In fact, (<&]_, $ 2 ) are second class constraints 

, f"(a 2 )a 3 Ni 

f*^> K 4ffer' < 3 - 17) 

Thus, the constraints can be solved provided we use Dirac brackets. In particular, the 
Hamiltonian takes a simple form 



1 3 2nip\ 

2 l 3 a? 



H = H=-a 3 1 a 2 n 1 --^, (3.18) 



where 

a 2 



Moreover, Dirac brackets for the variables ai,ni,7Ti,pi remain canonical. Therefore, 
eqs. (I3.18P and ( 13. 191) give the complete Hamiltonian description. We have checked ex- 
plicitly that it leads to correct equations of motion. In the case under consideration our 
formalism, when compared with Ostrogradski version, seems to be more complicated. 
However, it has an advantage that the curvature R is one of basic variables. 

4 Field theory 

Our formalism has a straightforward generalization to the field theory case. For simplic- 
ity, we consider only the Lagrangian densities depending on first and second derivatives. 
Such a density can be written in the form 

C = C($,d k $,d k d l $,$,d k $,ti>). (4.1) 

Again, we put $ = $i and select a function F = F(<&i, $i, $ 2 ) obeying 

d 2 F 
<9$i<9<l> 2 
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in the case of multicomponent field the relevant matrix should be nonsingular. We 
define, as previously, the function 

X = x(®i,dk$i,dkd l $ 1 ,$ l ,dkQi,$2), (4.3) 

as the (locally unique by virtue of (14.2J1 ) solution to the equation 



dx a$i 

Finally, the new Lagrangian density reads 



0. (4.4) 



£ = £($!, dkQudkdtQi, $1, &$i, X (- ■ •)) + g $ $1 

+ aF(<& 1 ,$ 11 <& a ) . a + aF($ 1 ,$ 1 ,$ a ) ' 

<9$ 2 <9<t>i 

It is now straightforward to check that the Lagrange equations 

^~ 4 a(w +4a <a(a^ (4 ' 6) 

yield the initial equation for the original variable $ = $1; as in the Section ITTTI $ = 
x{- ■ •)• One can now perform the Legendre transformation. The canonical momenta 
read 

ni (x) = —. , 1=1 d 3 xC. (4.7) 



Equations f)4.7|) can be solve (due to (14.21) ) with respect to $j: 

$! = $!($!, $ 2 , n 2 ) (4.8) 

<i> 2 = $ 2 ($ 1; 3 fe $ 1; ^^$!, 9^5 ro $!, $ 2 , 4$ 2 , 4^$ 2 , n 1; n 2 , <9 fc n 2 , 0^n 2 ). (4.9) 

if is defined in a standard way 

# = J (PxilL^xfaix) + n. 2 (x)& 2 {x)) - L, (4.10) 

and leads to the correct canonical equations of motions. 

Appendix 

A Extension to the case of arbitrary high deriva- 
tives 

Here we generalize the approach proposed in Section CD to the case of Lagrangians 
containing time derivatives of arbitrary order [21] . We restrict ourselves to the case of 
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one degree of freedom. We start with the Lagrangian depending on time derivatives 
up to some even order 

L = L(q,q,q\...,q^), (A.l) 

which is assumed to be nonsingular in Ostrogradski sense, - ,„ L , 2 7^ 0. Define new 

Qq(2n) 

variables 



^g (2j " 2) , 



,n + 1, 



Qi = 1 



(2i-l) 



l,...,n, 



so that 



L = L(qi,qi,q 3 ,q 2 ,...,q n ,q n ,q n+1 ). 

Let further F be any function of the following variables 

F — F{Ql, 9lj • • • ) Qn, Qn, q n +l, q n +2, • • • , ?2n) 

obeying 



+ 



Sg n+ i % 



0. 



and 



det 



d 2 F 
dqidqj 



7^0, n>2, 



«+2<i<2n 
l<j<n— 1 



(for n = 1 only ( 1A.5I) remains). 

Finally, we define a new Lagrangian 



„ r v^ /<9F. dF \ ^ 8F . 



Let us have a look on Lagrange equations 

dC _ d (d£\ 
dqi dt \dqi) 

Using ([Oil . (El), fTA~7l) and flA~8l) one finds 



0, i = l,...,2ra. 



«9 2 F 



^1 ^^ 



— (gfc+i - &) = « = n + 1, ...,2n. 



(A.2) 

(A.3) 
(A.4) 
(A.5) 
(A.6) 



(A.7) 



(A.8) 



(A.9) 



d 2 F 
dqidqj 



Consider the matrix 

n+l<i<2n 
l<j<n 

&f = for i = n + 2, . . . , 2n while jf^- 

oqidqn dq n +lOq n 



entering the LHS of eq. (IA.9|) . By virtue of (1A.5I) . 



a 2 L 



9ii + 



2 — 7^ due to Ostrogradski 



nonsingularity condition. Therefore, the first column of our matrix has only one non 
vanishing element. This, together with the condition ( 1A.6J) implies that it is invertible. 
Therefore, eq. ( 1A.9J) gives 

Qk+i = <7fe, k = l,...,n. (A. 10) 



Let us now consider (IA.8I) for 1 < i < n. We find 

dL d_ fdL\ _dF_ __ ^ /sf\ 

% dt\dqij dqi-i dt 2 \% / 

where, by definition, |f^ = 0. By combining these equations and using ( 1A.5I) and ( 1A.10I) 
we arrive finally at the initial Lagrange equation. We conclude that, contrary to the 
case of Ostrogradski Lagrangian, our modified Lagrangian leads to proper equation 
of motion. Let us now consider the Hamiltonian formalism. Again, the Legendre 
transformation can be immediately performed; neither additional Lagrange multipliers 
nor constraints analysis are necessary. In fact, let us define the canonical momenta in 
a standard way 

Pr = ^, (A.12) 



so that 



dqi 



dF 
Pi = — , i = n+l, ...,2n (A.13) 



Pi 



dL y, f d 2 F . d 2 F 

d 4i {^i \dqidq k k dqidq k 



2n Q2 



j=n+l ■> 



Due to nonsingularity of q q- . 1<i<2 „ eqs- (IA.13p can be solved for gi, q 2 , . . . , q n 

l<j<n 

Qi = fi(Ql,---,Q2n,Pn+l,---,P2n), 1=1,..., 71. (A.15) 

Now, eqs.l l A. 14]) are linear with respect to qi, i = n + 1, . . . , 2n and can be easily solved. 
Finally, the Hamiltonian is calculated according to the standard prescription. 

In order to compare the present formalism with the Ostrogradski approach let us 
note that they must be related by a canonical transformation. To see this we define 
new (Ostrogradski) variables qk,Pk, 1 < k < 2n: 

Qx-i = Qu i = l,-.,n, (A. 16) 

q 2 i = fi(qi,...,q2n,Pn+l,---,P2n), 1 = 1,..., 71, (A.17) 

dF 

P~2i-1 =Pi~ o— (gi,/l(- •• ),..., Qn, fn{-.. ),?„+!,..., g 2 n), 1= l,...,n, (A.18) 

dF 
P2i = --wriqijii- ■■),■■ .,?„,/„(...), fti+iv)fei)i z = l,...,n. (A.19) 

It is easily seen that the above transformation is a canonical one, i.e. the Poisson 
brackets are invariant. It is not hard to find the relevant generating function 

$(qi, ■ ■ ■ , q2n,Pl, q2,P3, &, ■ ■ ■ ,P2n-l, q2n) (A.20) 

n 

= ^2 QkP2k-i + F(qi, q 2 , q 2 , <?4, ■ ■ ■ , q n , hn, q n +i, ■■■, fen)- 

fe=i 
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Let us now consider the case of Lagrangian depending on time derivatives up to 
some odd order 

L = L(q,q,q,...,q^). 

Again, we define 



q t = q {2l - 2 \ i = l,...,n + l, 
ftsgC"- 1 ), z = l,...,n + l, 



so that 



L = L(q 1 ,q 1 ,q 2 ,q 2 ,...,q n+1 ,q n+1 ). 
Now, we select a function F, 

F — F(q±, q±, q 2 , $2, ■ ■ • , q n , q n , Qn+i, ■■■, ?2n+i) 
subject to the single condition 

8 2 F 



det 



dqidq k 



n+2<i<2n + l 
Kk<n 



7^0, 



and define the Lagrangian 



Consider the Lagrange equations (1A.8|) . First, we have 



E 



d 2 F 



k ^ dqidq k 
and, by virtue of (IA.26J) 



(?fc+i ~Qk) = Q, i = n + 2,...,2n + l, 



qk+i = qk, k = l,...,n. 

The remaining equations read 

dL d fdL\ dF d 2 ( dF 



TTT- + 



0, i = 1, . . . ,n + 1, 



c% dt \dqij dqi-i dt 2 \d<ji 
with fg = 0, gg^ = 0. Combining flA~29|) and flA~30|) one gets 



2n+l 






fc=0 



cft fc V^ (fc v 



0. 



Let us note that no condition of the form (1A.5I) is here necessary. 



(A.21) 



(A.22) 
(A.23) 

(A.24) 

(A.25) 

(A.26) 



" f dF dF \ 2n+1 dF 

£ ^ + £(%' l+ ^ +i ) + ,? +1 ^ (A - 27) 



(A.28) 



(A.29) 



(A.30) 



(A.31) 
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Also in the odd case the present formalism is related to that of Ostrogradski by a 
canonical transformation. Indeed, the canonical momenta read 



Pi 

Pi 



dF 

dL 

dqi 



+ £ 



n + 2,...,2n+l, 



d 2 F . d 2 F 

-Qk + r, . r, . Qk+1 



(A.32) 



2n+l 



+ £ 



j=n+l 



fc=l 

d 2 F 

dqidq. 



dqidq k dqidq k 



-Qj, 



l,...,n+ 1; 



(A.33) 



by virtue of (1A.26J) one can solve eqs. AA.32J) for qi,---,q n - The remaining n + 1 
equations (IA.33J) are used to compute the velocities q n+ i, . . . , (fen+i- m f ac t, using eqs. 
(1A.24I) - (1A.26J) as well as the Ostrogradski nonsingularity condition „.^ L 7^ one 
easily finds 

dPi ^ ^ 0. (A.34) 



det 



dqi 



l<i<n+l 
n+l<j<2n+l 



In particular 

Qi = fi(qi,---,q2n+l,Pn+2, ■■■,P2n+l), I = 1, . . . , fl. (A.35) 

Now, one can define the canonical transformation to Ostrogradski variables 



?M-1 = % 



l,...,n + l, 



(A.36) 

(A.37) 



q2i = fi(qi, ■ ■ ■ , q2n+l,Pn+2, ■ ■ ■ ,P2n+l), i = I,..., 71, 

dF 
P2i-1 =Pi~ — {qi, fi(. ..),... ,q n , fn(- ■ -),qn+i, ■■■ ,q2n+i), i = l,...,n+l, (A.38 



P2i 



dF 



dqi 

(qiJi 



Qn, fn{- • ■)) Qn+1, ■ ■ ■ j ?2n+l) 



l,...,n. 



(A.39) 



The relevant generating function reads 



$(gi, g 2 , • • • , q2n+l,Pl, q2,P3, fa, ■ ■ ■ , C[2n,P2n+l) 

n+1 

= 2_^P2k~iqk + F(qi, g 2 , • • • , On, <?2n, q n +li ■ ■ • , <?2n+l) 



(A.40) 



fe=i 



Summarizing, we have found a modified Lagrangian and Hamiltonian formulations 
of higher-derivative theories. They are equivalent to the Ostrogradski formalism in 
the sense that on the Hamiltonian level they are related to the latter by a canonical 
transformation. However, the advantage of the approach presented is that the Legendre 
transformation can be performed in a straightforward way. 

Acknowledgement Thanks are due to Prof. P. Kosihski for interesting dis- 
cussions. We are greatful to the unknown referee for very useful remarks. 
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